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Topological Excitations in Spinor Bose-Einstein Condensates 
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We investigate the properties of skyrmion in the ferromagnetic state of spin-1 Bose-Einstein 
condensates by means of the mean-field theory and show that the size of skyrmion is fixed to 
the order of the healing length. It is shown that the interaction between two skyrmions with 
oppositely rotating spin textures is attractive when their separation is large, following a unique 
power-law behavior with a power of -7/2. 
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Recent experimental successEr in obtaining Bose- 
Einstein condensation (BEC) of 23 Na atoms in an optical 
trap has opened up new directions in the study of BEC. 
Although confining atoms in a magnetic traps has been 
a successful method, it aligns the orientation of spin due 
to the existance of strong magnetic field. For this reason, 
spin of the trapped atoms cannot be regarded as a de- 
gree of freedom. In contrast, optical trap is based on the 
optical dipole force. Therefore atoms are confined in all 
hyperfine states. This offers the possibility of studying 
condensates with a spin degree of freedom, the spinor 
condensates. 

A variety of novel phenomena has been predicted for 
spinor condensates, such as changes in ground state 
structures with interaction parameters, propagation of 
spin waves and internal vortex structure.BBu) In par- 
ticular, the possibility of creating vortex states with- 
out core, or jSkyrmions, in spinor condensates has been 
pointed out.otP Skyrmions, which do not have an or- 
dinary vortex core due to the spin degree of freedom, 
offer us many new physical phenomena beyond those 
presented by other vortex states. One interesting fea- 
ture is the reduction of kinetic energy associated with 
the rotation by transferring this energy to the spin and 
eliminating energetically unfavorable normal core. This 
reduces the energy of skyrmion to a value lower than that 
of an ordinary singular vortex. 

In this paper we investigate the properties of 
skyrmions in spin-1 Bose condensates. First we discuss 
the basic properties of a single skyrmion, then we exam- 
ine the interaction between two skyrmions with opposite 
phases. 

We shall consider bosons with hyperfine spin F = 1. 
This includes alkali atoms with nuclear spin 7 = 3/2 
such as 23 Na, 39 K, and 87 Rb. The order parameter of 
the Bose condensate is characterized by 

(</>«(r)) = * Q (r) (1) 

= V^WCa(r) (a = 1,0,-1) 

where a indicates three hyperfine components, ip a ( r ) is 



the field annihilation operator of hyperfine state |l,a), 
n(r) = I'l'c^r)! 2 , ( a (r) is a normalized spinor i.e. 

OL 

Q ■ C = 1- This determines the average local spin 

(F) = C(r)F af3 (p(r). (2) 

HoB showed that in the case of spin-1 bosons consider- 
ation of two ground states is necessary since the effective 
interaction between two spins can be either ferromag- 
netic or antiferromagnetic. The ground states of the for- 
mer is minimized for (F) 2 — 1, and the latter, which 
is known for the polar state, is minimized for (F) = 0. 
The properties of vortices are different between these two 
states, and skyrmions are possible only in the ferromag- 
netic state. In this paper, we shall discuss their shape, 
size and interactions in the ferromagnetic state. 

The energy of the condensate is 

+ (C/-/i)n+^ 2 n 2 |, (3) 

where /x is the chemical potential, 52 = 47r7i 2 a2/M, ct2 
is the s-wave scattering length for two colliding atoms 
with total spin 2, and U is the trapping potential. 
All spinors in the ferromagnetic state are related each 
other by gauge transformation e l6 and spin rotations 
R(a,P,r) — e ~ lF * a e~ lF y l3 e~ lFz ' r , whece (a, f), r) are the 
Euler angles. The spinor £ is given byai> 

((r)=e* e R(a,[3,T) foj 

/i e -i(.-y+<*)(i + coa p)\ 

= " T72 e " l7sin ^ . W 
\ i e -^- Q )(l - cos/3) / 

where 7 = t— 9. The spin texture and supcrfluid velocity 
of the ferromagnetic state are given by 

(F) = (sin/3cosa)& + (sin (3 sin a) y + (cos f3)z, (5) 
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(V7 + cos/3Va). 



(6) 



Equation (5) shows that a, (3 characterizes the spin ori- 
entation. Moreover eq.(6) shows that spatial variation of 
spin in general leads to superflows. Note that the spin 
texture with superfluid velocity takes the same form as 
the angular momentum texture with superfluid velocity 
in 3 He-A.B> 

To illustrate the skyrmion, we use cyrindrical co- 
ordinates (r, <j), z) and assume that the spin texture and 
superfluid velocity of the condensate are both cylindri- 
cally symmetric. Taking a = 0, 7 = —<j), [3 — (3(r), 
where (3{r) is a function increasing from (3 — at r = 0. 
The spin texture and superfluid velocity field are 

(F) = (sin f3 cos 4>)x + (sin/3sin0)y + (cos (3)z. (7) 
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cos f3)e 4 



(8) 



Equation (8) shows that the superfluid velocity vanishes 
continuously as r — > CL and energetically unfavourable 
normal core is avoided.© For (3{r) — n/2, the superfluid 
velocity of the skyrmion reduces to an ordinary singular 
vortex. 

To obtain a precise order parameter of the skyrmion, 
we have to determine n(r) and (3{r) that minimize the 
energy. To investigate only properties of the skyrmion 
we neglect the trapping potential. From variational con- 
ditions 5K/5n — 0, 8K/5P = 0, we can derive a pair of 
coupled equations for n(r) and (3{r) 
2 
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Fig. 1. The variational function /3(r) when the spatial depen- 
dence of density is neglected. /3(r) is a monotonically increasing 
function obeying /3(0) = and /3(oo) = 27T. 



where R is the characteristic dimension of the condensate. 




da; 2 



2dAd/3 
/ dx J dx 



r(2-cos/3)sin/3 = 0, (10) 



where x — r/£ and / = \J n/n. Here we define ^ = 
(H 2 /2Mg2n) 1 / 2 as the healing length of the condensate 
and n = [if gi as the density of the condensate without 
flow. 

In view of the absence of the core, we shall neglect 
the spatial dependence of density for the moment, and 
replace it with an average n. In this case, the solution of 
eq.(10) is 

We plot this function in Fig. 1. Here we put /3'(0) = k. 
X = £/k expresses the typical size of the skyrmion. 

The energy of skyrmion per unit length along the z- 
axis is 



£s = 



(f) 4 +((f) 4 + 64)tan-(I(f) 2 ) 
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(R > A) 



(12) 
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Fig. 2. The energy of the skyrmion in case of neglecting spatial 
dependence of the density. For R 3> A, it approaches the con- 
stant (1 + it/2) 4Trh 2 n/M. 



It is clear from eq.(12) that the skyrmion has less energy 
than that of a singular vortex which is logarithmically 
dependent on the size of system due to the existence of 
singularity at the core. The R dependence of the energy 
of the skyrmion is shown in Fig. 2. 

Figure 2 shows that the energy of the skyrmion is 
minimum when A is infinite. This result is in fact not 
unreasonable physically. Our previous assumption has 
neglected the spatial dependence of density, and only 
considerd spatial dependence of £. Since £ tends to vary 
as slow as possible in order to decrease the kinetic en- 
ergy, the size of the skyrmion becomes infinitely large. 
However as the size of the skyrmion increase, the re- 
gion in which n(r) deviate from n becomes large. This 
increases the interaction energy. Therefore considering 
spatial dependence of density, the balance between the 
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kinetic energy and the interaction energy of the conden- 
sate determine the typical size of the skyrmion. 

We can assume without loss of generality that devi- 
ation of n(r) from the average density n is small i.e. 
n(r) = n(l — e(r)), e(r) <C 1 for r — ► oo. Linearizing 
eq.(9) with e(r) and using eq.(ll), we obtain the follow- 
ing equation 



dx 2 



-—-2e 

x dx 



= -2 



8k 3 x 2 + 4 3 fc 



fc 4 a 



, 4 3 



1 



(13) 



The solution of the above equation has the following 
asymtotic behavior 



e(x) 



de 



-V23 



C 2 Io(VW+2x) + -4 



2k 2 



(x> 1) 
(*« 1), 



(14) 



where Iq is a modified Bcssel function. Equation (14) 
shows that the variation of the density is spreaded over 
a distance £. Therefore the size of skyrmion is of the 
order of the healing length of condensate i.e. A ~ £. The 
skyrmion is thus stabilized by the fact that as it becomes 
large the gradient of spin texture on one hand decreases 
the kinetic energy and the variation of density increases 
the interaction energy on the other. 

Let us consider now two skyrmions with spin textures 
oppositely rotating 

i(l + cos/3(r)) 
&(r,0=| ^sin/3(r) 



Cs(r,0) 



±e 2i< ^(l - cos/3(r)) 

f(l + cos/3(r)) 
__2_ e ^ sin/ j( r ) 

i e 2^(i_ CO s/3(r)) 

4r/A 



(15) 



where (3{r) — 2tan~ 1 ( 8 _ 4 ^//^ ). The minus sign in front 
of the second component of £ g represents opposite phase. 
The texture in the presence of two skyrmions is schemat- 
ically represented in Fig. 3. 

We shall determine two-skyrmion energy as a function 
of their separation i? - Here we neglect the spatial de- 
pendence of density. The difference between the energy 
at a given separation i? and the energy of two isolated 
skyrmions per unit length is 

h 2 r 

£two(#o) = J d 2 r(VCtwo) t (VCtwo) - 2^.(16) 

For R 3> A we approximately express the order param- 
eter in the presence of two skyrmions as a sum of each 
order parameter. 



Ctwo iXi 



Cs(r_,<Ai) + Cs(r+,</>2). 



(17) 



Our result is plotted in Fig. 3. The curve has a mini- 
mum at i?o = ^0 where ro is approximately 5A. It can be 
seen that for Rq < ro the interaction between skyrmions 
is repulsive, while Rq > ro the interaction is attractive. 
As two skyrmions are brought together the region where 
spin texture varies gradually overlap, thereby changes 
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Fig. 3. (a) Schematic representation of two skyrmions. Bold ar- 
rows represent magnetic moment. The direction of circulation 
of each skyrmion is represented by a circular arrow, (b) Inter- 
section of the spin texture, (c) View from the z axis. 
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Fig. 4. The energy of two skyrmions as a function of the distance 
between the center of two skyrmions. 



the kinetic energy. Since texture tends to vary smoothly, 
skyrmions with large separation tends to come closer to 
each other, thus interaction becomes attractive. In con- 
trast, if they are too closer together, the region where 
velocity increases become large and the energy of two 
skyrmion increases. Thus with sufficiently close separa- 
tion the interaction energy is repulsive. 

We have investigated numerically the attractive inter- 
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Fig. 5. Logarithmic plot of energy of two skyrmions. The solid 
line, with a slope of -7/2 is the best fit to the data. 



action between skyrmions. Figure 5 shows logarithmic 
plot of energy of two skyrmions. It shows that attractive 
interaction between skyrmions is approximately propor- 
tional to Rq 7 ^ 2 ■ 

In conclusion, we have investigated the energy and the 
size of a skyrmion in ferromagnetic state of spin-1 Bose- 
Einstcin condensates and have shown that the size is of 
the order of the healing length. We also have shown that 
the skyrmion has lower excitation energy than that of 
a singular vortex. Therefore skyrmions can be excited 
more easily than singular vortices in scalar superfluids. 
This means that skyrmions are important in determin- 
ing thermodynamic properties of the condensate at low 
temperatures. It also means a reduction of the critical 
supcrfiuid velocity 

We have shown that the interaction between two 
skyrmions with opposite phases is attractive when their 
separation is larger than tq and repulsive when their sep- 
aration is smaller than r . Physical origin and the power 
of attractive force between skyrmions is still unknown, 
and interactions between skyrmions with other textures 
is also an interesting subject to invetigate. 
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